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Abstract 

We study the semi-classical trace formula at a critical energy level 
for an /i-pseudo-differential operator on K n whose principal symbol has a 
totally degenerate critical point for that energy. This problem is studied 
for a large time behavior and under the hypothesis that the principal 
symbol of the operator has a local extremum at the critical point. 

1 Introduction. 

The semi-classical trace formula for a self-adjoint /i-pseudo-differential operator 
P/i, or more generally /i-admissible (see |17|V studies the asymptotic behavior, 
as h goes to 0, of the spectral function: 

|Aj(/i)--E|<£ ' 

where the Xj(h) are the eigenvalues of Ph and where we suppose that the spec- 
trum is discrete in [E — e, E + e], some sufficient conditions for this are given 
below. If po is the principal symbol of Ph we recall that an energy E is regular 
when VpoO, ^ on the energy surface E B = {(x, £) e T*W n / p (x, £) = E} 
and critical when it is not regular. 

A classical property is the existence of a link between the asymptotics of (JIJ , 
as h tends to 0, and the closed trajectories of the Hamiltonian flow of po on 

lim ~?(E, h, ip) ^ {(t, i,()eIxl] E / $ t (x, = (x, 0}, 

where $t = exp(i7J po ) and iJ Po = d^po.d x — d x po.d^. A non-exhaustive list of 
references concerning this subject is Gutzwiller 10 , Balian and Bloch [I] for the 
physic literature and for a mathematical point of view Brummelhuis and Uribe 
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US], Paul and Uribe ^5], and more recently Combescure, Ralston and Robert 
|7], Petkov and Popov ^B], Charbonnel and Popov [H]. 

The case of a non-degenerate critical energy of the principal symbol po(x, £), 
that is such that the critical-set £(po) = G T*M™ / dp (x,£) = 0} is a 

compact C°° manifold with a Hessian cPpo transversely non-degenerate along 
this manifold, has been investigated first by Brummelhuis, Paul and Uribe in 
2 . They treated this question for quite general operators but for some "small 
times" , that is it was assumed that is the only period of the linearized flow in 
supp(c£>) when it is small. Later, Khuat-Duy in and ^2] nas obtained the 
contributions of the non-zero periods of the linearized flow with the assumption 
that supp(</?) is compact, but for Schrodinger operators with symbol £ 2 + V(x) 
and a non-degenerate potential V. Our contribution to this subject was to 
compute the contributions of the non-zero periods of the linearized flow for 
some more general operators, always with (p of compact support and under 
some geometrical assumptions on the flow (see 0] or 0). 

Basically, the asymptotics of Q can be expressed in terms of oscillatory 
integrals whose phases are related to the classical dynamics of po on . For 
(£o;£o) a critical point of po, it is well known that the relation: 

Ker(^$ f (x o ,£o)-Id)^{0}, (2) 

leads to the study of degenerate oscillatory integrals. What is new here is that 
we examine the case of a totally degenerate energy, that is such that the Hessian 
matrix at our critical point is zero. Hence, the linearized flow for such a critical 
point satisfies d x ^&t(xo, Co) = Id, for all i e R and the oscillatory integrals we 
have to consider are totally degenerate. 

The results obtained here are global in time, that is we only assume that 
supp(<y3) is compact. The core of the proof lies in establishing suitable normal 
forms for our phase functions and in a generalization of the stationary phase 
formula for these normal forms. 

2 Hypotheses and main result. 

Let Ph — Op™(p(x,£,h)) be a /i-pseudodiffcrcntial operator in the class of the 
^.-admissible operators with symbol p(x, £, h) ~ ^Pj( x , £), i- e - there exist 
sequences (p^j 6 E^(T*M") and (R N (h)) N such that: 

P h = h j pj{x, hD x ) + h N R N (h), \/N e N, 

j<N 

where Rn(H) is a bounded family of operators on L 2 (M. n ), for h < ho, and: 

S™(T*R") = {a : T*R" -> C, sup |<9£df a(x, £)| < C a ,pm(x, f), Va,f3 e N"}, 

where m is a tempered weight on T*R™. For a detailed exposition on h- 
admissible operators we refer to the book of Robert Let us notepo(a ; ,C) the 
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principal symbol of Ph, Pi(x,£) the sub-principal symbol and <£> t = exp(tH po ) : 
T*W l T*R n , the Hamiltonian flow of H po = d ( p -d x - d x p Q .d 6 . 

If E c is a critical value of po, we study the asymptotics of the spectral 
distribution 

1{E C M= E vi^^l (3) 

X j (h)e[E c -e,E c +e] 

under the hypotheses (-Hi) to (#4) given below. 

(H\) The symbol of Ph is real. 
There exists Eq > such that p^" 1 ([i? c — Eq,E c — £q\) is compact. 

Then, by Theorem 3.13 of JTj the spectrum a(Ph) H [E c — e, E c + e] is discrete 
and consists in a sequence Xi(h) < X2Q1) < ■•• < Aj(/i) of eigenvalues of finite 
multiplicities, if e < Eq and h is small enough. To simplify notations we write 
z = (x, £) for any point of the phase space. 

(H2) On the energy surface X.e c = Pq ^^({Ec}), po has a unique critical point 
z q = (xoi^o) and near zq : 

N 

p (z) = E c + ^Pj(z) + 0(\\(z - z )ir +1 ), k > 2, 

j=k 

where the functions pj are homogeneous of degree j in z — zq. 
(H 3 ) We have 6 Cft°(R). 

With (H2) the oscillatory-integrals we will have to consider are totally degen- 
erate. Hence, they cannot be treated by the classical stationary phase method. 
To solve this problem we impose the following condition on the symbol: 

(H4) The critical point Zo is a local extremum of pq. 

Remark 1 (H4) implies that the first non-zero homogeneous component pk is 
even and is positive or negative definite and also that zq is isolated on T,e c ■ 

Since we are interested in the contribution to the trace formula of the fixed 
point zq, to understand the new phenomenon, it suffices to study: 

Jz (E c ,h) = -^Tr / e^ l f(t)^ w (x,hD x )exp(-Up h )e(P h )dt. (4) 

Z7T J II 

R 

Here G is a function of localization near the critical energy surface £e c and 
ip € C§° (T*M. n ) has an appropriate support near zq. Rigorous justifications are 
given in section 3 for the introduction of <d(Ph) and in section 4 for i/j w (x, hD x ). 
Then, the new contribution to the trace formula is given by: 

Theorem 2 Under hypotheses (Hi) to (H4) we obtain: 

N 

lzo (E c , h) ~ h^- n (£2A j>k (<p)hi + O(h^)), ash^O, 
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where the Aj^ are some distributions and the leading coefficient is given by: 

Ao,feM = ^(^(i+Pi(^)),t? ri )^ r J \p k (9)\-^d6, (5) 

with t ZQ = max(t, 0) if Zq is a minimum and t ZQ = max(— t, 0) for a maximum. 
Remark 3 One can derive a full asymptotic expansion as shows Lemma H5I 

3 Oscillatory representation. 

Let be <p e S(R) with 6 C °°(M). We recall that: 

1 (E c ,h)= v( Xjih) ~ Ec ), I s = [E c -e,E c + e], 

with Pq 1 (I Eo ) compact in T*R™ the spectrum of Ph is discrete in I e for £ < £o 
and h small enough. Now, we localize near the critical energy E c with a cut-off 
function 6 £ Cg°(}E c -e,E c + e[), such that 9 = 1 near E c and < 9 < 1 on 
R. The associated decomposition is: 

7 (£ C) ft)= 7l (£ c ,ft)+7 2 (£; c ,ft), (6) 

yi(E c ,h)= £ (i-e)(M%(^), 
l2 {E cl h)= Q(XAh)M Xj{h) - Ec ). 

The asymptotic behavior of "fi(E c , h) is given by: 
Lemma 4 7i(J5 c , ft) = e>(ft°°) as ft -> 0. 

Proo/. Since if E <S(R), Vfc G N, 3C fe such that |sc fc <p(a:)| < C fe on E. By Theorem 
3.13 of |I7| the number of eigenvalues N(h) lying in I E n supp(l — 9) is of order 
0(h~ n ), for ft small enough. This gives the estimate: 

hi(E c ,h)\<N(h)C k \ Xj{h) - Ec \- k . 

But on the support of (1 — 9) we have \Xj(h) — E c \ > So > 0. This leads to: 

\li(E c ,h)\ < N(h)C N eo k h k < c N h k - n . 
Since the property is true for all k £ N, this ends the proof. ■ 



with: 
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Consequently, for the study of j(E c ,h) modulo 0{h°°), we have only to 
consider the quantity 72(^0 h). By inversion of the Fourier transform we have: 

R 

Since the trace of the left hand-side is exactly 72(^0 h), we obtain: 

y 2 (E c ,h) = -^Tr / e^( / 3(i)exp(-^P h )e(P /l )di, (7) 
Z7T j n, 

R 

and with Lemma 0] this gives: 

j(E C) h) = ^Tr / e ^^Wexp(-^p i )e(P h )^ + 0(^ 00 ). 

R 

Remark 5 Another interest of this formulation is that, under the geometrical 
condition to have a "clean" flow, j(E c ,h) can be expressed, at the first order, 
as the composition of two Fourier integral-operators. 

Let be Uh(t) = exp(— ji-Ph) the evolution operator. For each integer L we 
can approximate Uh(t)Q(Ph), modulo 0(h L ), by a Fourier integral-operator, 
or FIO, depending on a parameter h. To give a precise formulation of this 
approximation we recall briefly the principal notions on FIO. Let be N G N, 
Lp € C°°(R" x R N ) and a € C£°(R™ xR N ). An oscillatory-integral with phase 
ip and amplitude a is 

I(ae^) = (2tt/i)-t J a(x,e)e^ v(x - e) d9. (8) 

To attain our objectives, it suffices to consider amplitudes a with compact sup- 
port. We suppose, as usually, that tp — <p(x, 9) is a non-degenerate phase func- 
tion, i.e. 

d(d 8l <p) A ... A d(de N ip) ^ on <t(cp) = {(x, 6) e R" x R N / d 9 (p{x, 6) = 0}. 

This implies that £(<p) is a sub-manifold of class C°° of R™ x M. N and that: 

. f <%) -> T*(R"), 
** : \ foflJ.-fod^M)), 

is an immersion. In this situation one say that ip parameterizes the Lagrangian 
manifold A v = i v ((L(ip)). Conversely, if A C T*R ra is a Lagrangian sub-manifold 
we can always find locally some non-degenerate phase functions parameterizing 
A, see e.g. |Hj- Now if ipi and ip2, <Pj & C°°(R™ x R^), are two non-degenerate 
phase functions parameterizing locally the same Lagrangian manifold, i.e. A^ n 
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U = A V2 n U, for an open U C T*R", then it is well-known that there exists a 
constantc G R such that for all ai G Cg°(R n xR Nl ) with i Vl (supp(oi)nC(v?i)) C 
[/, small enough, there exists a sequence a-i j G C^°(R n x R" 2 ), j G N, such 
that for all L G N : 

I{ ai e^) = e^ c J2 h 3 I(a 2 ^ei V2 ) + h L r L (h), 

with rj,(/i) uniformly bounded in L 2 (R"), since a\ G C^°, and c comes from 
— S V2 = c on A Vl n {/ = A V2 fl U. This recalls that the fundamental object 
associated to an oscillatory-integral is the Lagrangian manifold parameterized 
by the phase function. The next definition follows Hormander. 

Definition 6 Let be A C T*R" a Lagrangian sub-manifold of class C°° . The 
class I(R n , A) of oscillatory functions associated to A is given by: 

u h G J(R n ,A) &u h = Y J I{av{h)e^), 

where the sum is locally finite and the functions a v (h) are of the form: 

a v (h)= h j cij, dj G C°°(R n x R Nv ). 

-d v <j<j v 

With d v , J v and N v in N and where ip v is a non degenerate phase function 
parameterizing A n i Vv (supp(a„(/i))), where supp(a v (h)) — (J supp(a t , J ). 

3 

Now let be A C T*(M™ x R ( ) a C°° Lagrangian sub-manifold of T*(R n xl 1 ). 

Definition 7 The family of operators (F h ) : L 2 (R") L 2 (R l ), < h < I, is 
a h-FIO, associated to A, if there exists two families and R\^\ WeN, 

such that: 

i) Each F^ has an integral kernel in I(R n x R l , A). 

ii) For each N the operator is bounded on L 2 and there exists Cm > 
such that \ \R^ \\c{L 2 (R n ),L 2 (R 1 )) < Cn, uniformly in h, < h < 1. 

in) For each N : F h = F^ N) + h N R { h N) . 

After these definitions, we recall the theorem that gives the approximation 
of Uh{t)<d(Ph)- Let A be the Lagrangian manifold associated to the flow of pn- 

A = {(t, r, x, £, y, r,) G T*R x T*R n x T*M. n : r = p (x, £), (x, = r?)}. 

Theorem 8 TTie operator Uh(t)@(Ph) is an h-FIO associated to A, there exists 
C^efcW w ^ integral kernel in 7(R 2 ™ +1 ,A) and R\f\t) bounded, with a L 2 - 
norm uniformly bounded for < h < 1 and t in a compact subset of R, such 
that U h (t)e(P h ) = U ( ^(t) + h N R[ N \t). 
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We refer to Duistermaat 8 for a proof of this theorem. 

Remark 9 By a theorem of Helffer and Robert, see e.g. ^7], Theorem 3.11 
and Remark 3.14, 0(-F\) is an h- admissible operator with a symbol of compact 
support in p ( ^ 1 (/ e ). This allows us to consider only oscillatory- integrals with 
compact support. 

For our goal the following corollary is crucial. 

Corollary 10 Let be 9i G Cq°(M.) such that<di = 1 on supp(6) and supp(Oi) C 
I e> then VN e N: 

Tr(0(P^(^=^)) = ^Tr J ^e^U^Q^dt + 0(h N ). 

R 

Proof. By cyclicity of the trace, for all N € N we have: 
Tr(6(P,M^^)) = ±Tr J ^^{U^it) + h N R { h N) {^{P^dt. 

R 

But Qi(Ph), and a fortiori Qi(Pf l )R < j^\t) by the ideal property, is a class-trace 
operator (see e.g. ^7J, section 2.5), with norm: 

We^R^m^ < \\R { h N) (t)\\ 1101(^)11^ < c^mu, 

since it is true for all t € supp((^) the corollary is proven. ■ 

The geometrical approach associated to the Hormander class I(R 2n+1 , A) 
gives a great degree of freedom in the choice of the phase function, this will be 
exploited below. In fact if (x ,£o) 6 A and if ip = <p(x,6) <E C°°(R" x R N ) 
parameterizes A in a neighborhood U, small enough, of (xo,£o) then for each 
u h e I(W\A) and x e Cg°(T*R n ), supp(x) C U, there exists a sequence of 
amplitudes dj — aj(x, 9) <E C^°(R n x R N ) such that for all integer L: 

X w (x,hD x )u h = hU{ aj e^) + 0{h L ). (9) 

-d<j<L 

We will use this remark with the following result of Hormander (see ^JJ, tome 
4, proposition 25.3.3). Let (T,r,xo,t;o,yo,—r/o) G Afl ow , r/o ^ 0, then near this 
point there exists, after perhaps a change of local coordinates in y near y$, a 
function S(t, x, rf) such that: 

<P(t,x,y,r)) = S(t,x,rj) - (y,ri) , (10) 

parameterizes Afl ow . In particular this implies that: 

{(t,d t S(t,x,ri),x,d x S(t,x,ri),d v S(t,x,ri),-r))} C A flow , 
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and that the function S is a generating function of the flow, i.e. 

^t{d v S{t,x,rf),rf) = {x,d x S(t,x, V )). (11) 

Moreover, S satisfies the Hamilton- Jacobi equation: 

f d t S(t,x,rj) + p (x,d x S(t,x,rj)) = 0, 
\ S(0,x,£) = (x,O- 

Now, we apply this result with (xq, £o) = (yo, rjo), our unique fixed point of the 
flow on the energy surface Ee,.. If £o — we can replace the operator P^ by 
e%( x >^> p h e~T;( x 'iu with £i ^ 0. This will not change the spectrum since this 
new operator has the symbol p(x, £ — £i, h) and the critical point is now (xo, £i) 
with £i ^ 0. Consequently, the localized trace J2(E C , h), defined by Eq.JfjJ), can 
be written for all N € N and modulo 0(h N ) as: 

j 2 (E c ,h) = J2( 2 ^ h )~ d+j I e^ s ^ x ^-^ +tE ^ aj {t,x,^{t)dtdxd^. 

j<N Rx „ 2 „ 

(12) 

To obtain the right power — d of h occurring in Eq. lQ we apply results of 
Duistermaat jS] (following here Hormander for the FIO, see tome 4, for 
example) concerning the order. An /i-pseudodifferential operator obtained by 
Weyl quantization: 

(Jhrfc)"* J a(^±y,t)ei(*-y<$dt, 

R N 

is of order w.r.t. 1/h. Now, since the order of Uh(t)Q(Ph) is — |, we find that 



^ w (x,hD x )U h (t)e(P h ) ~ 5^(27rft)- n+ J / a J (i,x,y,r;)et( s ( t - a; ''')-^'''»dy, 

i<iv R J „ 

(13) 

Multiplying Eq. (|13fl by <p(t)e"K tEc and passing to the trace we find l|12l) with 
d = n, where we write again dj(t, x, if) for Oj(t, x, x, if). 

To each element uh of 7(K n , A) we can associate a principal symbol a pi i nc (uh) , 
where S is a function on A such that £dx — dS on A. If Uh = I(ae^' p ), then we 
have S = S v = <p o z^ 1 and cr pr i nc (-u/ l ) is a section of |A| 5 <g> M(A), where M(A) 

is the Maslov vector-bundle of A and \A\i the bundle of half-densities on A. If 
pi is the sub-principal symbol of Ph, the half-density of the propagator Uh(t) is 
easily expressed in the global coordinates (t, y, if) on Aflow via: 

t 

exp(i J p x {® s {y,-rf))ds)\dtdydi)fi. (14) 
o 

This expression is related to the resolution of the first transport equation for 
the propagator. For a proof we refer to Duistermaat and Hormander [H]. 
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4 Classical dynamics near the critical point. 

A critical point of the phase function of Eg. (|12l leads to the equations: 

E c = -d t S(t,x,0, , , , 
x = d s S{t,x,0, <sJ .77 
t = £s\t\x,& 1 **(».0 = (*.0, 

where the right hand side defines a closed trajectory of the flow inside £,e c . 
Since we are interested in the contribution of the critical point, we choose a 
function ip G C^°(T*R™), with tp = 1 near z , hence: 

j 2 (E c ,h) = i-Tr f e i! ^0(t)iP w {x,hD x )eiq>(-^tP h )Q(P h )dt 

R 

+ ^Tr j £(t)(l - ip w (x,hD x ))exp(-^-tP h )G(P h )dt. 

R 

The contribution on supp(l — ip) is non-singular and can be treated by the 
regular trace formula. If supp(V0 is small enough only points (t, x, £) = (t,Zo), 
for t £ supp(<p), will give a contribution, since zq is isolated on Y>e c - 

Now, we restrict our attention to the contribution of the critical point. Until 
further notice, the derivatives d will be taken with respect to initial conditions. 
Since zq is totally degenerate, we obtain: 

d^t(za) = exp(O) = Id, V*. (15) 

With (H2), the next homogeneous components of the flow are given by: 

cP $t(*o) = 0, Vt, Vj G {2, .., k - 2}. (16) 

To obtain the next non-zero terms of the Taylor expansion of the flow, we will 
use the following technical result: 

Lemma 11 Let be zq an equilibrium of the C°° vector field X and $ t the flow 
of X . Then for all m £ N*, there exists a polynomial map P m , vector valued 
and of degree at most m, such that: 

t 

m—lff* I ~ \( ~ra — 1\ 



d m <£ t (z a )(z m ) = d$ t (z ) / d<£^ s (z )P m (d<£ s (z )(z),...,d m - L <S> s (z a )(z m - L ))d S . 

(17) 







Proof. We note x l G (R ra ) z the image of x under the diagonal mapping, with the 
same convention for any vector. If /, g are two C°° applications, by the "Faa di 
Bruno formula" , we have: 

d m (fog)(x ){x m ) = C( V )cTf(g(x ))(d^g(xo)(x^),...,d n ^(x )(x n n), 
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where T(m) is the set of partitions of the integer m and where it is assumed 
that the partition p of m is given by m = n\ + ... + n r and C(p) are some 
universal integers. Hence, for our fixed point zq we obtain: 

d m (Xo<P s )(z )(z m ) = C-(p)<rA-(z )(d ni *a(«o)(2 ni ),...,d nr *.(ab)(« T,r )). 

For Y — (Yi, Y m ), we can define: 

Pm<X)= C{p)d r X(z )(Y ni ,...,Y nr )-dX(z )(Y m ), (18) 

p£P(m) 

this leads to the differential equation, operator valued: 
^(d m $ s (zo))(z m ) = dX(z )d m $ s (z )(z m )+P m (d® s ^ 

With the initial condition d m $o(- 2 o) = 0, the solution is given by Ea. (|17fl . ■ 

Since d<&t(zo) — Id, Vt, the first non-zero term of the Taylor expansion is: 

d fc " 1 $ t (zo)(z fc ~ 1 ) = J d^HpMiz^ds = td^H^zo)^- 1 ), (19) 
o 

where the last identity is obtained using that d fc - 1 J Br Po (zo) = d k - 1 H Ph (z ). 
Moreover, with crpo(zo) = 0, for the next term Lemma ITT1 gives: 

t 

d k $ t (z )(z k ) = J d k H po (z )(z k )ds = td k H Pk+1 (z )(z k ). 
o 



Remark 12 For the derivatives dP<&t(zo)i with j > k there is two different kind 
of terms, namely: 

t 

J diH Po (z )(z j )ds = tdiH PQ (z )(zi) = 0(\\z\\i), 
o 

and terms involving powers of t. For example, we have: 




+2 



= L rf ,+2- fejffpo(zo)(zJ+ i- fc;dfc -i iJpfc(zo)(zfc -i )) ^ 

This term is simultaneously 0(t 2 ) and C(||^| P) near (0, zq). A similar result 
holds for other terms, which are 0(t d ) for d > 2, by an easy recurrence. 
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Lemma 13 Near z$, here supposed to be to simplify, we have: 

S(t,x,£) - (x,S) + tE c = -t(p k (x,Z)+R k+ i(x,Z)+tG k+1 (t t x,Q), (20) 

where R k+1 (x,t) = O(\\(x,0\\ k+1 ) and G k+1 (t,x,0 = O(\\(x,0\\ k+1 )> uni- 
formly for t in a compact subset o/M. 

Proof. By Taylor and under hypothesis (H2) we obtain: 

* t (x,t) = (*,£) + —L^d*- 1 *,^)^*- 1 ) + 0(\\z\\ k ). (21) 
Now, we search our local generating function as: 

N 

s(t, x , = -tE c + (x, o + Y, s *(*> x > + ociifo oir +1 ), 

i=3 

where the functions Sj are time dependant and homogeneous of degree j w.r.t. 
(x,£). With the implicit relation: $t(d£S(t,x,£),£) = (x, d x S(t, x, £)) and using 
Eq.J2U we have: 

S(t,x,0 = -tE c + (x,0 + S k (t,x,0 + O(\\(x,0\\ k+1 ), 

and comparing terms of degree k — 1 gives: 

jvs k (t,x,n) = -—l—d^StCoxfoO*- 1 ), 

where J is the matrix of the usual symplectic form. By homogeneity and with 
Ea. lfT^jl we obtain that: 

s k (t,x,0 = l- l ({x,z),tJd!'- 1 H Ph (x,z) k - 1 ) = -tp k (x,0- 

It remains now to treat the remainder. Since S(0, x,£) = {x, £), we have: 

S(t,x,Z)-(x,Z)=tF(t,x,0, 

with F smooth in a neighborhood of (x, £) = 0. Now, the Hamilton-Jacobi 
equation imposes that F(0,x,£) = — Pq(x,£) and we obtain: 

R k+1 (x,0= P0 (x,Z)-E c -p k (x,0 = O(\\(x,0\\ k+1 ). 

Finally, the time dependant remainder can be written: 

S(t, x, - S(0, x, - td t S(0, x, = 0{t 2 ), 

since by construction this term is of order 0(||(:e,£)||' c+1 ) we get the result. ■ 



11 



5 Normal forms of the phase function. 

Since the contribution we study is local, we can work with some coordinates 
and identify locally T*R™ with R 2 ™ near the critical point. We define: 

*(t, z) = *(t, x, = S(t, x, - (x, + tE c , z = {x, e M 2 ". (22) 



Lemma 14 J/ P/, satisfies conditions (H2) and (H4) then, in a neighborhood 
of z = Zq, there exists local coordinates x such that $f(t,z) ~ X0X1 */ z o is a 
maximum and ^(t,z) ~ — X0X1 if z o is minimum. 

Proof. We can here assume that zq is the origin and we use polar coordinates 
z = (r, 6), 6 € 8 2 " _1 (R). With Lemma IT3l near the critical point we have: 

¥(M)~¥(t,rd) = -tr k {p k (6) +rR k+1 {9) +tG k+1 (t,r9)), 

where p k (9) is the restriction of p k on S 2n_1 . We define new coordinates: 

(X0,X2, ...,X2n)(t,r,8) = (t,dl, ...,6» 2 „-l), 

Xi (t, r, 0) = r |p fc (0) + rR k+1 (9) + tG k+1 (t, r0)\* . 

In these coordinates the phase becomes — X0X1 f° r a minimum and X0X1 f° r a 
maximum. Now, we observe that: 

^(t,O,0) = \p k (9)\K Vi, 
hence, the corresponding Jacobian satisfies the relation: 
\J X \(t,0,0) = \pk(9)\i ^0, Vi, 

and this defines a local system of coordinates near z$ for all t. ■ 
This change of coordinates leads to: 



: ^(xo,Xi)dXodxi 



with a new amplitude A defined by: 

Mxo,Xi)= [x*(a(t,r6)r 2n - 1 \Jx\)d X 2...d X 2n. (23) 



Remark 15 Since Xi{t,r,6) = r|j)*(0) + rR k+1 (9) + tG fc+ i(i,r6>)|fc, our new 
amplitude satisfies A( X o,Xi) = 0(Xi n )i near Xi = 0- This will play a major 
role since distributions of Lemma Il6l below are supported in {^i = 0}. 
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We end this section with a lemma on asymptotics of oscillatory integrals. 

Lemma 16 There exists a sequence {cj)j of distributions, whose support is 
contained in the set {xi — 0}, such that for all function a £ C^°(IR+ xR): 



( / e iX ^a( Xo ,Xi)dxo)d X i -^A^c^a), (24) 



(i 



asymptotically for A — > oo, with 

Cj = ~(F(x 1± ^ ± )(xo)®5 { 3) (xi)), x- =m&x(-x,0). 

K J ■ 

Proof. We note (t,r) for (xo>Xi) ancl we define g(r,r) = ^ r t (a(i, r))(r), where 
J- t is the partial Fourier transform with respect to t. Then, we obtain: 



( J e lXtrk a(t,r)dt)dr = J g(-Xr k ,r)dr = \~i J g(-r k , -^)dr. 

OK 

A Taylor expansion with respect to r for g(r, r) at the origin gives: 



Afc z=o 



where Rn+i (t, A) is integrable with respect to r, with L 1 norm uniformly bounded 
in A. By a new change of variable we obtain: 

oc 

if, i r 1 v — 1 i+i f d a , , , , (+i-fc , N+i , 

X ~~ k J 9(~r k ,jr)dr= ^]i X ~^ / (r,0)|r|-^dr + 0(A~^). 

n 1 — 



If we introduce X- = max(— a;, 0) and T(x_ k )(r) the result follows. ■ 

Remark 17 A similar result holds for the phase — XoXi if we change xo into 
— Xo and we then have simply to replace terms X- by x + = max(x, 0). 

6 Proof of the main result. 

We can assume that zq is a maximum. Lemma [TBI shows that the first non-zero 
coefficient is obtained for I = 2n — 1 (sec Remark H~5ll and is given by: 

l (2n- 1)! (^(^)^ S o^ 1) )^(Xo,Xi)) = \ yV(^)(Xo)A(xo,0)d X o 
where, by construction, the new amplitude is: 

i(xo,Xi) = / x*(a(t,r0)\Jx\\Pk(O)+rR k+1 (9)+tG k+1 (t,r0)\-^)dx2-dx2n. 
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Now, we use an oscillatory representation of A(xo,0) vi a: 
Mxo,0) = ^ / A( XQ ,Xi)e lzXl dzd X i. 



2tt . 

We return to the initial coordinates and, since xo — i> we obtain: 

Mt,0) = ^ f a(t,r9)\p k (9) + rR k+1 (d) +tG k+1 (t,r6)\-^e tz ^ t ^ 9 '>dzdrd9. 
2tt J 

Inserting the definition of %i, we use the change of variable: 
y = z\p k {6) + rR k+ i(0) + tG k+1 (t, r6)\i. 
Since in r = we have Gk+i(t, 0) = 0, by integration in (y, r) we obtain: 

A(xo,0)= / a(xa,0)\H(0)\-^d9. (25) 



zn — K 

The distribution T(x_ k )(xo) of Lemma ITr)l acts on this function via: 

(H^)(xo),Mxo,0)) = (^(aC^),a(.,0)) J \p k (6)\-^d6. 

g2n-l 

Hence, the top-order contribution to the trace formula is given by: 

7(^c,/ l ) = ^7^ TT ^(^r^) ) a(.,0)) y Ipfc^l-Tde + O^-^-*). 

With a(t,0) = 0{t) exp(itpi(zo)) , see (|X4f> . by Fourier inversion we have: 
^(^(x^),a(.,0)) = J <p(t + Pl (z Q ))t^ dt. 

R 

Finally, if zq is a local minimum we must simply replace i_ by t + = max(i, 0) 
and this complete the proof of Theorem [21 I 
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